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Abstract 



Recently, the E821 experiment at the Brookhaven National Laboratory 
announced their latest result of their muon g — 2 measurement which is about 
2.6-0" away from the standard model prediction. Taking this result seriously, 
we examine the possibility to explain this discrepancy by the supersymmetric 
contribution. Our analysis is performed in the framework of the unconstrained 
supersymmetric standard model which has free seven parameters relevant to 
muon g — 2. We found that, in the case of large tan/3, sparticle masses are 
allowed to be large in the region where the SUSY contribution to the muon 
g — 2\s large enough, and hence the conventional SUSY search may fail even 
at the LHC. On the contrary, to explain the discrepancy in the case of small 
tan/3, we found that (i) sleptons and SU{2)l gauginos should be light, and 
(ii) negative search for the Higgs boson severely constrains the model in the 
framework of the mSUGRA and gauge-mediated model. 
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The latest result of the E821 experiment at Brookhaven National Laboratory on the 
muon anomalous magnetic moment (7^ — 2 

a^(E821) = (g^ - 2)/2 = 11659202(14)(6) x 10"^° (1) 

indicates a possible confrontation with the standard model of particle physics. The deviation 
between the experimental value and the standard model prediction is 

a^(E821) - a^(SM) = 43(16) x 10~^°, (2) 

or 2.6-cr deviation [||]. In order to fill up the discrepancy, a new physics beyond the standard 
model is called for. The apparent deviation is comparable to or even larger than the contribu- 
tion from the standard model electroweak sector computed as a^(SMEW) = 15.1(0.4) x 10~^° 
0. This suggests that the energy scale of the new physics should be very close to the elec- 
troweak scale and/or it should have some enhancement mechanism to give a large contribu- 
tion to a^. 

Among other things, supersymmetry (SUSY) is the most promising candidate for such 
a new physics. To solve the naturalness problem in the Higgs sector of the standard model, 
the superparticle masses should lie below the TeV scale. The SUSY contribution to has 
been investigated in the literatures (see [0HTO[| and references therein) .Q Generically it is 



sizable for superparticles weighing less than 1 TeV. And as we will explain shortly, it is 
enhanced for large tan /3 region, where tan (3 is the ratio of the vacuum expectation values 
(VEVs) of the two Higgs bosons in the SUSY standard model. 

The purpose of this paper is to reexamine the contribution to in the framework of the 
unconstrained minimal supersymmetric standard model (MSSM) in the light of the recently 
reported experimental data. The unconstrained MSSM has more than 100 SUSY breaking 
parameters and one usually impose some relations among the model parameters: otherwise it 
would be very difficult to trace all dependence of the parameters for some specific processes. 
The SUSY contribution to a^, however, depends only on seven MSSM parameters as we 
will list below, and thus we can leave them as free parameters to analyze their dependence. 
Another important point that makes the model independent analysis possible is that, unlike 
flavor changing neutral current (FCNC) processes, we do not rely on particular mechanisms 
to suppress SUSY FCNC to compute the SUSY contribution to a^. Conclusions we will 
draw are therefore very general. 

The apparent deviation from the standard model prediction implies a non-vanishing 
SUSY contribution. We will identify the parameter region of the MSSM which is capable to 
account for the discrepancy (0) at 2-a. Exphcitly we require the SUSY contribution to lie 
in the following range, 

11 X 10-^° < a^(SUSY) < 75 x 10"^°. (3) 

The SUSY contribution to decreases as the superparticles become heavier. Since a non- 
vanishing SUSY contribution is needed, we will obtain an upper bound on the mass scale 



^See Ref. [11,12| for scenarios of large extra dimensions. 
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of the superparticles. For large tan (3 region, the enhancement mechanism works and @ is 
easily satisfied even for relatively large superparticle masses. We will find that the Wino as 
heavy as 1 TeV can be compatible with it. Another point we pay a particular attention is 
the case of low tan (3. We will show that for tan /? ^ 3 the SUSY contribution can be large 
enough to explain the deviation, without confronting the present bounds on the superparticle 
masses obtained by negative searches of superparticles at collider experiments. 

At one-loop level, the SUSY contribution to stems from chargino-sneutrino loops as 
well as from neutralino-smuon loops. Formulae of the SUSY loop contributions are given, 
for example, in Ref. ||^. In the (unconstrained) MSSM, the parameters involved with this 
process are the soft SUSY breaking mass parameters for the right-handed and left-handed 
smuons denoted by m^^ and m^^, respectively, the trilinear scalar coupling for muon A^, 
the U{1)y and SU{2)i gaugino mass parameters Mi and M2, the higgsino mass parameter 

and tan [3 which is the ratio of the VEVs of two Higgs bosons. Throughout this paper, 
we take the parameters to be real, and do not include possible CP phases which are known 
to be generically small due to stringent limits on electric dipole moments. (See, however, 
Ref. [|13|.) For most of the analysis we impose the GUT relation for the gaugino masses 
Ml : M2 ~ 1 : 2. We checked that our following results are almost unchanged even if the 
GUT relation is relaxed. In the MSSM analysis, we do not consider any particular SUSY 
breaking scenario, and thus we do not introduce specific relations among the remaining 
six parameters mentioned above. For comparison, we will briefiy discuss the case of the 
minimal supergravity scenario (mSUGRA) later on. Furthermore we do not impose the 
lightest superparticle (LSP) in the MSSM sector to be neutral. A charged LSP would be 
ruled out if it were stable, but there are many ways out, including R-parity violation and 
light gravitino LSP. 

For generic SUSY mass parameters, it is known that the chargino-sneutrino diagram gives 
a dominant contribution to a^(SUSY). Then the relevant parameters to calculate a^(SUSY) 
are m^^, M2, /i and tan/?. Thus it will be reasonable to fix the other mass parameters for 
some specific values for the moment. Here we take = and assume the GUT relation 
to the gaugino masses.0 The sign of the SUSY contribution to is directly correlated with 
the sign of M2/i in most of the parameter regions. It is positive (negative) for M2/i > 
(M2/i < 0). Thus, the result of the E821 experiment given in Eq. (|l]) suggests M2IJ, > 0, and 
hence we consider the positive sign case in the following. 

Another important point is that the SUSY contribution to is enhanced for large tan f3 
because, in the dominant diagrams, muon chirality is fiipped by the muon Yukawa coupling 

oc 1/ cos P ~ tan/5, not by the muon mass itself Thus, for the large tan/? case, 

the SUSY contribution can be large enough to explain the discrepancy even with relatively 
heavy superparticles. On the contrary, when tan f3 is small, Wino and slepton masses should 
be light to make the SUSY contribution to large enough. 

In the framework of the unconstrained MSSM, we calculate an upper bound on the 
lighter smuon mass m^i as a function of the Wino mass parameter M2. The result is shown 
in Fig. |I[ In deriving the upper bound, we vary m^^ and m^^ , and derive the largest possible 



^We checked that the value of a^(SUSY) is insensitive to the value of unless is extremely 
large. 
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value of mp_i which can reahze the 2-cr bound (0). In Fig. |T], we take tan/3 = 3, 5, 10, 30, 
and we fixed fi = 500 GeV as a typical example. We checked that the lighter smuon is 
mostly the left-handed one which consists of the SU{2)l doublet with the sneutrino as it is 
involved with the dominant chargino-sneutrino loop. We find from this figure that quite a 
large parameter region is in accord with the consideration of the recent data of a^. Thus we 
conclude that the SUSY is naturally able to explain the apparent discrepancy observed at 
the experiment. 

Fig. |l| also shows that the constraint on the lighter smuon mass m^i is much stronger 
than that on the Wino mass M2. For instance, for tan/5 = 10, the smuon must be lighter 
than about 400 GeV while the Wino as heavy as 1 TeV is allowed at 2-cr. 

Let us now closely look at the large tan/3 case. In this case, the SUSY contribution 
is enhanced as tan/? increases. Thus one can expect that even heavy superparticles can 
be compatible with the lower bound a^(SUSY) = 11 x 10~^°. In Fig. ||, we show a plot 
of the 2-0" upper bounds on m^i as a function of the /x parameter. Here we take several 
values of tan j3 and M2 = 1 TeV. Even with such a large Wino mass, we find that the 2-a 
constraint can be satisfied with tan/3 as small as 10 when the slepton mass is lighter than 
about 200 GeV. Applying the GUT relation of the gaugino masses to the gluino mass as 
well, the Wino mass M2 = 1 TeV corresponds to the gluino mass of about 3 TeV. Namely 
the accuracy of the present data still allows the possibility of such a heavy gluino in 
the framework of the unconstrained MSSM. The Large Hadron Collider (LHC) experiment 
will not be able to reach such a heavy gluino, and thus SUSY searches at the LHC would 
require unconventional approaches. A future analysis of with more statistics will further 
reduce the error in the measurement, which may constrain the superparticle masses within 
the reach of the LHC. 

Let us turn to the case of low tan /3. In this case the enhancement mechanism operated 
at large tan/3 is not effective and thus one may expect that the SUSY contribution is 
rather small. The case of tan /3 = 3 in Fig. 1 shows, however, that the SUSY contribution 
can explain the deviation at 2-a level when the SUSY mass parameters are close to their 
experimental bounds. 

Next we would like to discuss the mass of the lightest scalar Higgs boson rrih- It is known 
that small tan (3 tends to give relatively light mh because the tree-level contribution to m\ 
is roughly given as m| cos^ 2/3. To survive the Higgs mass bound obtained at LEP 200, 



which approaches 113.5 GeV as the pseudo-scalar Higgs mass increases [^, large radiative 



correction from top and stop loops is required [|r^. Thus the stop masses should be large 
enough. Roughly speaking, the mass bound is satisfied for small tan /3 ~ 3 if the stop masses 
exceed about 1 TeV. The relation between the slepton masses which are constrained by the 
analysis and the squark masses is highly model dependent. 

To give more quantitative arguments on the light Higgs mass, we consider the case of 
the mSUGRA. In Figs. ^ and ||, we plotted contours of O-a, l-cx, and 2-a preferred values of 



^Existence of light superparticles may affect the fit to the electroweak precision data, which was 
considered in Ref. [|lj] in the parameter region where a^(SUSY) is sizable. According to Ref. |14|, 
there is some region of light higgsino dominant LSP case in which the inclusion of SUSY particles 
gives a better fit to the electroweak precision data than the standard model alone. 
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a^(SUSY) on mo vs. M2 plane, where mo is the universal scalar mass at the GUT scale and 
Ml/2 is the universal gaugino mass at the same scale. (The Wino mass at the electroweak 
scale is related to M1/2 as M2 ~ 0.83Mi/2.) In this framework, we also calculated the 
lightest Higgs mass m/i, and plotted the contours of constant m^ on the same figures. In 
the figures the trilinear scalar coupling is taken to be zero. As one can see, for the case of 
tan/5 = 5, the constraint uih > 113.5 GeV severely constrains the parameter region which 
gives preferred value of a^(SUSY). Taking account of the Higgs mass constraint, however, 
the 2-0" constraint (^ can be realized with tan/3 ^ 5. It is interesting to see that the region 
with small universal scalar mass is favored. Although the region with mo < 100 GeV is 
not shown in the figures, we checked that 2-cr constraint and the Higgs mass constraint can 
be simultaneously satisfied in the limit of mo (i.e., with the no-scale type boundary 
condition). For smaller tan/3, a^(SUSY) cannot be large enough to explain the discrepancy 
at the 2-cr level. For large tan /3, on the contrary, the Higgs mass can be easily large enough 
in the region with sufficient a^(SUSY). 

Here, let us comment on the allowed parameter region in the mSUGRA case. We find 
that the allowed parameter region is rather tighten in the mSUGRA scenario. To illustrate 
this point, let us consider the case of tan/? = 5. In the unconstrained MSSM, the 2-a bound 
constrains M2 ^ 600 GeV. On the other hand, in the mSUGRA case, larger value of M1/2 
results in larger slepton mass through renormalization group effect. As a result, the 2-a 
bound in the mSUGRA case is M1/2 ^ 300 GeV which corresponds to M2 ^ 250 GeV.0 
The difference will be crucial when one considers the discovery potential of superparticles 
at hadron colliders. 

To make a comparison, we also consider a case of the gauge-mediated SUSY breaking 
(GMSB) scenario whose free parameters are the overall scale of the soft-breaking parameters 
Aqm, the messenger scale Mmess, tan/3, and the number of the 5 + 5 messenger multiplets 



iVmess [171. The overall scale is related to the Wino mass as 



M2 = ^^^Agm, (4) 

with g2 being the SU{2)l gauge coupling constant. In our analysis, we use M2 instead of 
Aqm using this relation. In Fig. ^, we show the region which is consistent with the 1-cr and 
2-cr constraint on M2 vs. tan/3 plane. Here, we take Mmess = 10^ GeV and A^mess = 1- In 
addition, we also plot the constant irih contour. We can see that a^(SUSY) can be within 
the 2-cr bound with tan/3 ^ 5 even with the Higgs mass constraint m^ > 113.5 GeV. We 
checked that this result is insensitive to the choices of Mmess and N^f^^s- 

Finally, let us briefly comment on the branching ratio of 6 — > 57. To this process, 
the SUSY contribution is dominated by the stop-chargino loops. It is well-known that it 
can interfere with the standard model and charged Higgs contribution constructively or 
destructively, depending on the relative sign of the parameter and the trilinear scalar 
coupling of stop At. Furthermore, in most cases the sign of At is essentially determined 



^We expect that a similar tight bound on M2 can be obtained in a wide class of models where the 
supersymmetry breaking effect is mediated at high energy scale, and thus the left-handed smuon 
mass suffers from the renormalization effect of SU{2)l gaugino. 
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by the sign of the gluino mass M3 via large renormahzation group effect. It is found that 
for Ma/i < the SUSY contribution is constructive, which is in fact disfavored because the 
branching ratio B{b 57) tends to be predicted too large compared to the experimental 
value. Thus the b s'y consideration favors the M3/X > case. On the other hand, the 
positive contribution to the is obtained when the sign of M2/i is positive, provided that 
the chargino loop dominates over the neutralino loops. Thus the case that the Wino and 
gluino have masses with the same sign is favored by the combined consideration of and 
B{b ^ sj). This is the case in many models of SUSY breaking. In particular, models with 
the GUT relation of the gaugino masses fall into this category, including the mSUGRA and 
the GMSB. On the other hand, the case with the opposite sign of the Wino and gluino masses 
is disfavored, which is in fact the case in a simple class of anomaly mediated supersymmetry 
breaking models |^. 

In this paper, we considered implications of the recently reported measurement at the 
E821 experiment to supersymmetric standard models. We made the analysis mainly based 
on the unconstrained MSSM. Since the SUSY contribution should be non-vanishing at 2-a 
level, the upper bounds on the superparticle masses were obtained. For large tan /3, the 
superparticles can be quite heavy, which may be escaped from the LHC gluino reach. We 
observed that the allowed region of the superparticles masses is significantly larger than the 
case of constrained models such as the mSUGRA model. On the contrary for smaller tan /3, 
the Wino mass as well as the lighter slepton mass should be light. In this case, the bound 
on the Higgs boson mass obtained at LEP200 gives stringent constraints. We illustrated 
this point in the mSUGRA model as well as the GMSB model, yielding tan /3 ^ 5 for both 
models. 

Given the upper bounds on the masses of the sleptons and the charginos/neutralinos, 
lepton flavor violation such as /i — 67 as well as fi-e conversion may be observed in near 



future experiments [jl8l. Another implication is to proton decay. As we discussed, a model 
dependent lower bound on tan/? is obtained by combining analyses of a^(SUSY) and nih- 
This tightens proton decay constraints in SUSY GUT models, though details are quite model 
dependent. Result of detailed study along this line will be presented elsewhere. 

To conclude, we should emphasize the importance of the further reduction of error of the 
measurement which is expected to be done in near future, as well as the further study 
of the uncertainty coming from the hadronic vacuum polarization. We hope that they may 
sharpen the confrontation with the standard model more clearly in near future and confirm 
the necessity of the physics beyond the standard model. 

While preparing the manuscript, we were aware of the papers: A. Czarnecki and 



W.J. Marciano, [hep-ph / (J 1 (J2 1 22|, L. Everett, G.L. Kane, S. Rigolin and L.-T. Wang, [hep- 



ph/0102145| , J.L. Feng and K.T. Matchev, |hep-ph/0102T46| , E.A. Bahz and R Gondolo, 
[hep-ph/0102147| , U. Chattopadhyay and P. Nath, |hep-ph/01021571 , which have some over- 
lap with our analyses. 
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FIG. 1. The upper bounds on the hghter smuon mass m^i. The horizontal hne is the Wino 
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and 30 from below. 



500 GeV and = 0. Here, tan/3 is taken to be 3, 5, 10, 



1000 



M2 = 1000 GeV 



> 

CD 
O 

w 

CO 

cc 
E 

ii 
o 

C 

o 
Si 

CD 
Q. 
Q. 
3 



100 



tan p =50 
tan p =30 
tan p =20 



Ntanp=10 



100 



1000 



II [GeV] 

FIG. 2. The 2-cr upper bounds on the lighter smuon mass m/ii. The horizontal line is the 
Higgsino mass parameter /x, and we take M2 = 1 TeV, = 0, and tan/3 = 10, 20, 30, 50. 
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FIG. 3. Contours of the constant a^(SUSY) on mo vs. M2 plane in the minimal supergravity 
model. The solid (dotted, dashed) line is for the center (±1-0", ±2-0") value of a^(SUSY). We take 

= and tan/3 = 5. We also plot the constant nih contours in the dash-dotted lines {mh = 110, 
115, and 120 GeV from below). 
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FIG. 5. Contours of the constant a^(SUSY) on M2 vs. tan/3 plane in the gauge-mediated 
SUSY breaking scenario. The solid (dotted, dashed) line is for the center (il-cr, ±2-(t) value of 
a^(SUSY). We take Mmcss = 10^ GeV and iVmcss = 1- We also plot the constant contours in 
the dash-dotted lines, {rrih = 110, 115, and 120 GeV from below). 
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